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Abstract. We consider a system in direct contact with a thermal reservoir and which, 
if left unperturbed, is well described by a memory-less equilibrium Langevin equation 
of the second order in the time coordinate. In such conditions, the strength of the 
noise fluctuations is set by the damping factor, in accordance with the Fluctuation and 
Dissipation theorem. We study the system when it is subject to a feedback mechanism, 
by modifying the Langevin equation accordingly. Memory terms now arise in the time 
evolution, which we study in a non-equilibrium steady state. Two types of feedback 
schemes are considered, one focusing on time shifts and one on phase shifts, and for 
both cases we evaluate the power spectrum of the system's fluctuations. Our analysis 
finds application in feedback cooled oscillators, such as the Gravitational Wave detector 
AURIGA. 
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1. Introduction 

In recent years, several studies (both experimental and theoretical) have analyzed 
cooling protocols of a signal generated by a thermal source [HO [3], HI [5]. This is achieved 
using a feedback apparatus attached to the measuring device, and should be intended 
as an effective cooling from the bath temperature T to the lower effective temperature 
T e ff. In essence, a fraction of the signal generated by the thermal source is fed back to 
counteract the input signal. This in turn has the global effect of decreasing the signal 
amplitude, as if the measured source has been further cooled. A possible application is 
that of quantum oscillators [6]- Typically, as of today, practical limits prevent one from 
achieving bath temperatures that are low enough for the quantum nature of the oscillator 
(e.g. a nanomechanical device) to be clearly detectable. Feedback cooling is used in 
the attempt to 'close the gap'. Another notable application is that of gravitational 
wave detection [3 E], such as in the AURIGA detector [1], where the feedback is used for 
technical reasons [9]. 

Our intent here is to develop a model of the behavior induced by a feedback 
apparatus on a resonant circuit whose current is driven by a stochastic voltage of thermal 
origin (see for instance Figured]). 

Some studies have already dealt with similar problems [TQl E], with one crucial 
difference. The feedback was represented by a delayed potential in a first order 
Langevin equation, which generates a characteristic oscillatory behavior. We stress 
that, differently, we focus on the case of a harmonic oscillator. Moreover, the feedback 
may act on any term of the Langevin equation. We also show the difference between 
applying the feedback as an ideal delay line in the time domain and a filter in the 
frequency domain, a distinction that becomes important in many actual realizations. 

In the following we shall: i) illustrate how the second order equilibrium Langevin 
equation should be modified in the presence of feedback terms of all orders in the time 
derivative; ii) consider the case of an RLC circuit in which the feedback driving is the 
time derivative of the current coupled to a memory term; iii) derive the power spectra 
of the current in stationary states when the feedback employs a delayed time-shift (the 
'digital protocol'), highlighting the fundamental differences that emerge in the case 
of high, moderate and low quality factors for the electronic oscillator; demonstrating, 
among others, that the resonance frequency in all cases depends on the feedback gain 
factor, and that the effect of changing the time shift is markedly different depending 
on the quality factor; iv) derive the above power spectra if the feedback employs a 
phase shift (the 'analog protocol'), a situation equivalent to that of AURIGA, choosing 
different cutoff frequencies and discussing similarities and differences with the digital 
protocol. 
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2. From equilibrium to feedback driving 

Let us review the case of a system in equilibrium with a thermal reservoir, which 
can be described by a one-dimensional equilibrium Langevin equation [12], 

q + 7^ + u%q = rj (1) 

where r] = r](t) is a Gaussian noise for which the following FDT relations must hold: 

(#)) = o 

(r,(t)r J (t'))=2 1 (q 2 )8(t-t f ) (2) 
2lkBT 8{t - t>) (3) 



m 
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The quantity (q ) is the square of the observable q(t), averaged over the statistical 
distribution of the noise (therefore a constant in stationary conditions). Equation 
(J2J) is a constraint that holds for any stationary states generated by the dynamics 
([T]), if r](t) is delta-correlated and satisfies the causality relations (r](t)q(t')) = and 
(r](t)q(t f )) = 0, for t > t'. Equation ([3]) involves thermodynamic equilibrium properties, 
presupposing that the quantity q represents one quadratic degree of freedom in an 
underlying Hamiltonian (m being the corresponding 'mass' proportionality factor), 
provided the canonical distribution at temperature T holds. 

Equation ([3]) encompasses a multitude of systems well described by a harmonic 
oscillator driven by a stochastic force. For the specific case of an RLC circuit coupled 
with a thermal reservoir, q(t) is the charge and I(t) = q(t) is the current. The circuital 
parameters R (resistance), L (inductance) and C (capacity) are such that: 7 = R/L, 
Uq = 1/LC, m — L. The thermal bath thus produces a stochastic voltage Vr(t) with 
correlations dictated by FDT, i.e. (Vr(t)VT(t')) = 2ksTR6(t — t'). This means that the 
resistance dissipates all the fluctuations induced by the thermal bath. 

We then connect the circuit to a feedback apparatus. The observable quantities at 
earlier times are continuously stored and fed back into the system at later times. We 
wait until new stationary states (if they exist) are reached. This corresponds to a new 
non-equilibrium situation generated by the constant driving influence of the external 
apparatus. Correspondingly, Equation (jHJ) is modified by the inclusion of the feedback 
terms. Here we deal with settings that lead to a modification of the dynamics that 
makes room for memory terms, while still preserving linearity, such as, 

K(t-t')q(t')dt / + [ \{t-t')q{t')dt'+ [ fi{t-t')q{t')dt' = pit) (4) 



If n(t) is kernel-invertible, we shall write, 

q(t)+ I 7(t - t')q(t')dt' + f p[t - t')q{t')dt' = N{t) (5) 

We assume that the stochastic processes pit) = p{r](t)} and N(t) = N{r](t)} are 
functionals of the same noise-generating source which operates at equilibrium, with 
no additional source of stochasticity. In other words, we consider a noiseless feedback. 
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I d (t) = F{I s (t)} 
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Figure 1. The scheme of a RLC driven by a thermal voltage Vt and connected to 
a feedback device which 'translates' the input current I s (t) into Id(t), a functional of 
I s {t). A particular realization of this scheme has been already implemented to analyze 
the output of AURIGA [HQS]. 



p(t) = p{r)(t)} and N(t) = N{r)(t)} are supposedly known from analyzing the feedback 
mechanism. Therefore, although Equation (jSJ) closely resembles well known Langevin 
equations with memory, it is not guaranteed that FDT-type relations (T2J [13] are directly 
applicable here. In typical situations, one arrives at formulas resembling Equation (JSJ) 
after a series of qualitative considerations and reasonable physical assumptions about 
how one should separate the many degrees of freedom into a deterministic evolution and 
a stochastic driving [Ti], At that point, causality conditions are applied, assuming the 
separation was done correctly. Here we face a rather different non-equilibrium situation, 
somewhat simpler to some extent. In fact, the functional expression for the cumulative 
stochastic process N(t) = N{rj(t)} is in principle given once the feedback mechanism is 
known explicitly. At the same time however, N(t) depends on the past history of r)(t), 
rendering the assumption of causality less justified. This issue is somewhat reminiscent 
of what has been discussed in other works [15] and will be clarified later. 



3. Concrete examples - types of feedback 

Let us now proceed with the concrete example of Figure [U The scheme works by 
first amplifying the input signal current I s (t). After an appropriate transformation, a 
fraction of the original current is then returned to the RLC. For such a transformation 
we consider two protocols, which we name the 'digital' and the 'analog' protocols. 

The equation describing the dynamics of the circuit is, 

(L - L m )q(t) + Rq{t) + -q(t) = V T {t) - L in q s (t) 

I s (t) = I(t) + I d (t) (6) 

where Vr(t) is the voltage generated by the thermal bath, and we have neglected the 
noise generated by the feedback apparatus. 

The current Id{t) is ceded back to the main circuit, i.e. it is the feedback output. 
It is a functional (here linear) of the current I s (t), which in turn is both the measured 
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signal and the input of the feedback. The dependence of Id(t) on I s (t) is determined 
by the feedback. It is helpful to express the currents in the frequency domain. Since 
we consider only stationary states, we can use Fourier transformation. The Fourier 
transform of f(t) is given by 



f(u) = / e^f(t)dt. 

J — oo 

In the digital protocol Id{t) is given by I s (t), scaled by a gain factor G < 1, and shifted 
in time by td- 

I d (t) = GI s (t-t d ). (7) 

The analog protocol is best described in the frequency domain, and we shall consider 
a low-pass filter acting on a current proportional to I s (t), 

AO 

W = FT— -M") (8) 

i L LUJ 

with Q the angular cut-off frequency and A the gain in the amplification line. 

Equations (J7|) and (jSJ) imply that the current I d (t) depends on the past history of 
the circuit, thus revealing that Equation ([6]) is indeed in the family of Equation (j4j). 

In either case, in Equation ([H]) one can eliminate I d (t), or its Fourier transform 
Id(u), in favor of I s (t) or I s (uj) respectively. 

The digital protocol is a recurrent scheme in theoretical studies. Literally, the 
feedback is treated as an ideal external 'field' which continuously stores and releases 
the dynamical entries, exploiting an ideal delay line. This model is better suited if 
one evaluates numerically the temporal evolution with a simple updating rule. In 
practice, however, it is not obvious whether such ideal systems actually arise in nature, 
and whether they are easy to be operationally implemented. That the feedback is 
rather obtained by a filter working in the frequency domain seems a valid alternative. 
Incidentally, the latter is precisely the way the gravitational wave detector AURIGA 
operates [I]. 

The question then arises as to when the two schemes are equivalent to each other, 
or in what way they differ and under what circumstances. As a matter of fact, in the 
case in which the currents are almost sinusoidal with angular frequency ujq because, 
e.g., of the oscillator's high quality factor, choosing td = ir/2uJo and Q <C uq leads to 
Id{t) — Gu q s (t) in both cases (where we took AQ = Gu ). It has already been shown 
that, to a good approximation, the feedback induces an extra dissipation of the order 
of L in Gu [HEB]. 

Consider now the two protocols. 

4. The digital protocol - I d (t) = GI s (t - t d ) 

Eliminating I s (t) from Equation ([6]), one has, 

Lq{t) - G(L - L in )q\t - t d ) + Rq{t) - RGq{t - t d ) 

+^q(t) - ^q(t - t d ) = V T (t) - GV T (t - t d ) (9) 
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This equation is of the type of Eq. (Ill) , and shows that the time evolution of I is 
determined both by the voltage acting directly at time t, and by its portion continuously 
fed back into the system. When G = 0, one recovers the equilibrium situation described 
by Equation Q, with rj(t) = V T {t)/L. 

We define x = L in /L, y — 1 — x, 7 = R/L, 1/LC = uj 2 . As t can be arbitrarily 
large in a stationary state, we can systematically eliminate q at earlier times in (JSj), to 
obtain a convergent series (if \Gy\ < 1) like, 



00 

x \ — \ 



N(t) = q( t )+ 1 q{t)+u 1 q{t)--^{Gy) *[ 7 <K*- HO +"o?(* (10) 

^ fc=i 

where iV(t) is given by 

00 

iV(t) = [V T (t) - - Y,(Gy) k V T (t - kt d )]-. (11) 
y k=i 

Often, the measured quantity is the signal current I s (t), rather than I(t). Then, 
eliminating I(t) in favor of I s (t) one has 

Lq s (t) - G(L - L in )q s {t - t d ) + Rq s (t) - RGq s (t - t d ) 

+^qs(t)-^q s (t-t d ) = V T (t) (12) 

and, analogously to Equation ffTUl) . 

N s (t) = q s (t) + iqs {t) + ulq s (t) (13) 

00 

x 



~ ~ Y,(Gy) k [iq s (t - kt d ) + LU 2 q s (t - kt d )\ 

y L 



k=l 



with 



N s (t) = j Y,(Gy) k V T {t-kt d ). (14) 

Equations (11011 lj) and fll3lll4p are equations of the type of Eq. ([3]). The unknown is 
the voltage describing the interaction of the circuit with the thermal reservoir. The 
distribution assigned to N(t) and N s (t) will induce the probability distributions of 
{/(*), <Z(*)} and of {I s (t),q s (t)}. 

Thus, to make progress we need some assumption about the statistical behavior of 
the voltage generated by the thermal bath. As customary, we assume that the driving 
due to the bath is a Gaussian delta-correlated process also when the feedback apparatus 
is switched on. We thus assume (Vr(t)Vr(t')) oc 8(t — t'). To fix the constant we assume 
that - as in equilibrium - all the fluctuations are dissipated by the resistance R, so that 
(V T (t)V T (t')) = 2k B TR5(t - t'). Alternatively, we shall write (rj(t)rj(t')) = 2v6(t - t'), 
with v = ksTR/L 2 = ksT'-f/'m. Formally, this implies that (fj (u) fj (uj')) = 4ttv5(uj~\-uj'). 
Still by Fourier transformation, one can calculate the power-spectra of the currents / 
and J, as, 



hoo 



Sj{u) = / e^(J(0)J(t))dt ; J = /, I s . (15) 



00 
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By transformation of both sides of Equation (TTU]) and summation over the index k, 
we can express conveniently the formal Fourier transform of the current I(t) (here 
a stochastic variable) as, 

ifa) = -icq^) = ^(1 _ Gye***) 

with 

D{u) = u 2 {l - Gye tu)td ) +iury(l- Ge iutd ) - cu 2 (l - Ge luJtd ). 

From ffl3l) . an equivalent expression holds for I s (u), with N s (uj) in place of N(uj). 
Equation (TlTj) yields 

A~r/ \ 1-Ge™* s 
N(u) = - - — — ■qfu) 

while ( Tl4l) yields 



N s {u) 



1 - Gye iu)td ' 

The above expressions also indicate that the transfer functions T(u>) = I(u)/fj(u) and 
T s (u) = I s (u})/fj(cJ) are readily available, with 

- Ge"*') 

™ = m 

and 



DM 



Then, recalling that J(t) = |^°°e %u)t J{uo)duj /2-k, via the hypothesis (77(0^)77(07') ) = 
47rw5(u; + a/), ffl5l) can be written as, 

„ , , 2v u; 2 (l + G 2 - Ge ia " d - Ge^) . . 

= D(uj)D(-u) (16) 

= n( 2 w \ (17) 
D{uj)D{—uj) 

Obviously, the alternative derivation, starting directly from equations ([§]) and (JT21) . leads 
to the same result. 

We focus on the expression ( fl7|) . since I s is typically the signal that is measured in 
an experiment [1]. The denominator reads explicitly, 

D(uo)D(-uo) = a 2 + G 2 (3 2 + (1 + G 2 ) 7 V (18) 
+2xG'ju 3 sin (ut d ) - 2G(a(3 + rfuj 2 ) cos (ut d ) 

with 

a = uj 2 — cJq, P = y u 2 — Uq (19) 
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Under the hypothesis that v = kBTR/L 2 as in equilibrium, Equation (1171) with 
the specification f|T8|) constitutes a quantitative estimate for the full spectrum of the 
current, if the low-noise feedback apparatus has an ideal delay line with gain G G (0, 1). 

Regarding the parameters of the circuit, we analyze different possibilities. First, 
we consider a high quality factor 7 <C ujq, i.e. R^C/L <C 1. We further assume that 
7 <C xGuo <C ujq. If the feedback operates as an effective cooling, or damping, it is 
appropriate to tune the delay time td by setting t d u ~ ir/2. uj q is the natural resonance 
frequency of the oscillator, defined as luq = 1/ y 'LC ' . 

In the proximity of the resonance frequency, t^uj — vr/2 implies that cos (ut^) ~ 
and sin (uitd) — 1. The term cubic in uj in Equation ( f!8l) is also negligibile. After 
neglecting the smallest terms, one verifies that the denominator in ffl7]) is approximately, 

D{uj)D{-u) ~ (1 + y 2 G 2 )u 4 - 2(1 + yG 2 )uo 2 uj 2 + (1 + G 2 )^ (20) 

which implies a power-spectrum of a Lorentzian form, 

2^ 2 

The resonance frequency of the circuit, i.e. the frequency that maximizes this power 
spectrum, is u r . At equilibrium, Z = v = k B TR/L 2 , n — 7 = R/L and u r = u = 
l/y/LC. With the feedback active, both the resonance frequency and the effective 
dissipation become dependent on the feedback gain. More specifically, 7 <^ xGu lj 
entails 

"t = ^k< (22) 



y 2 G 



.2 




1 + yG 2 



1 \ \l 1 I „.2^2 ~ I +y 2 G 2l U ° 



1 + y 2 G 2 

The expressions for the resonance frequency and for the damping factor can be further 
simplified when 1 — y — x — L in / L 1 and G C 1, leading to, 

xG 2 



U r + YVG 2 ) ^ ~ XGUJ ° (23) 

Thus, the resonant frequency explicitly varies with the feedback gain, albeit 
very slightly, a feature that may have contacts with other sit uat ions [TTj . The power 
spectrum of the measured current is essentially equivalent to that of a similar circuit 
at equilibrium, with resonance in the proximity of uq but with damping increased by 
a relative factor approximately xGujq/^ = L in Guo/R and temperature lowered by the 
same factor. The fact that the feedback can indeed act as an effective additional damping 
is well established [U [2j El H [5]. 

In Figure [2] we show the power spectra (fT7|) . using Q = luq/i = VL/(RVC) = 10 5 
and x = Li n /L = 1CT 2 . We choose four different values of G: 0.06, 0.15, 0.4 and 
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Figure 2. Effective cooling at very high quality factors. Vertical axis in logarithmic 
scale in base 10. The power spectrum of the signal current I s as a function of the 
normalized frequency lu/u>q. (I 2 ) e q is the mean squared current at equilibrium. Data 
are for ujq/j = Q = 10 5 ; x = Li n /L = 10~ 2 ; f</Wo = tt/2. From top to bottom; the 
curves are for G = 0.06; G = 0.15; G = 0.4; G = 0.75. The solid lines represent 
Equation (|17|) for the four settings. The (o) points are from (|2"Tj) and {jT^i . showing 
almost perfect overlap. The dashed lines represent the further approximations of Eqs. 
(|23|) . For the two top curves at G — 0.06 and G — 0.15 we only show one line because 
the three expressions are indistinguishable. Vertical line is set at uj = luq as a guide for 
the eye for the slight shift of the resonance uj r away from with varying G. 

0.75. We compare the full expression (I17j) with the curves deduced from the Lorentzian 
( 12Tj) . using respectively the set of values of Equation (1221) (the o points) and the set 
of Equation ( 1231) (the dashed lines). For G = 0.06 and G = 0.15 the three curves are 
almost identical, but their differences can be appreciated for G = 0.4 and G = 0.75. We 
conclude that (TPTj) approximates very well a Lorentzian with parameters (122|) . and to a 
lesser degree is also consistent with ( 1231) . 

In Figure [31 we show the peculiar behavior of the spectrum derived from Equations 
( 1T71) and ( [HD , when we set the delay time tdUJo = ir. We compare the two results with the 
spectrum obtained in the absence of feedback and with that obtained with i^o = 7r /2- 
Rather than the over-damping characteristic of t^uJo = 7r/2, = 7r entails a marked 
shift of the resonance frequency with respect to equilibrium. 

The whole picture changes if we start decreasing the value of Q. For example, let 
Q = 100. As shown in Figure HI in this case two features develop. A set of secondary 
peaks appears, showing that the spectrum ceases to be purely a Lorentzian and becomes 
an oscillatory function of its argument uj. The reduction of the spectrum amplitude 
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0.9995 0.99975 1 1.00025 1.0005 
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Figure 3. The power spectrum of the current I s at Q = 10 5 , for two different choices 
of the delay time, against the equilibrium spectrum. Coordinate axes are as in Figure 
[21 Represented are the equilibrium curve (the solid line), and two curves with active 
feedback with x = 10~ 2 ; G = 0.06. The o points show the shift t^coo = ir/2. The 
dashed line is for t^tuo = n. 

around the natural resonance frequency ujq is now accompanied by a slight amplification 
in correspondence of other frequencies. Furthermore, if we restrict our attention to 
the close proximity of the equilibrium resonance frequency (the inset of Figure HI), the 
setting tdU>o = n compares very differently than before to either the equilibrium case or 
the setting ts^JQ = tt/2. See for example Figure [3] and our previous discussion. These 
two features show that in these conditions, despite Q still being relatively high, the 
feedback drives the system away from its equilibrium configuration in a characteristic 
way, quite different from the previous case. 

Finally, Q — 1 sets a radically different scenario. Now the dissipation is so high 
that the feedback has a more dramatic impact on the frequency profile. As we see in 
Figure [5l changing the delay t^ makes the resonance frequency u r (that at which the 
spectrum has a maximum) to jump discontinuously from one value to another, for the 
same gain factor G. This is due to the fact that the closest secondary peak increases in 
amplitude with growing 

One can explore a variety of circuital parameters, assuming the general expression 
for the power-spectrum ([17)1 holds. For very high Q, we have identified the Lorentzian 
(T2~TT) and the parameters Z, u r and \i. This is a quasi equilibrium scenario, although 
radically different from it in a thermodynamic perspective (see [IS]). As the quality 
factor decreases, we have seen some special features develop. The feedback that operates 
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Figure 4. Effective cooling at moderate quality factors. Development of secondary 
peaks. Vertical axis in logarithm scale in base 10. The power spectrum of the 
signal current I s as a function of the normalized frequency lu/ujq. Data are for 
<V7 = Q= 10 2 ; x = L m /L = 1CT 2 ; G = 0.75. Thick solid line shows t d uj Q = ir/2; 
thin solid line is for t^uio = 7r; dashed line is the equilibrium case (absence of feedback). 
Inset: Zoom in around lvq. 



at a fixed frequency, as in the digital protocol, radically changes the frequency profile of 
the oscillator. The question is, how much of this richness is an artifact of the feedback 
operating with a fixed delay time? 

In the next section we shall try to answer to this question, by considering a different 
feedback mechanism. 

5. The analog protocol - /^(o;) = -^^-I s {u) 

As we have said, an alternative way to implement the feedback protocol is to let 
the amplified current AI s [t) be passed through a low-pass normalized filter of cutoff 
frequency H < w . In this case, Equation ^ is best managed directly in the frequency 
domain. We eliminate q(u) and I(u) in favor of q s (uj) and I s (uj), and in order to do 
that we introduce 



By Fourier transforming Equation ([6]) and using the usual notation, one finds, 





(24) 
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Figure 5. Effective cooling at low quality factors; the 'bimodar trend. The power 
spectrum as a function of the normalized frequency uj/loq- Data are for 7/wo = Q = 1; 
x = L in /L = 10~ 2 ; G = 0.75. Solid line shows td^o = 0.587T. +-dashcd curve is for 
tduio — 0.807T. The dashed line is the curve at equilibrium. The arrows indicate the 
trend of the two peaks with increasing td- The two peaks do not overlap during the 
'transition'. 



B(u) = uo 2 [(l-yA)Q-iuj} 
+ 270; [(1 - A)Q - iu] - - A)Q - iuo] 

Since I s (u) = —iuq s (uj), Equation (12^1) leads once again to an expression of the transfer 
function for the current, T s (u) = I s (uj)/fj(uj) = (oj 2 + iuQ)/ B(u), and the power 
spectrum now reads, 

SM = B(u)B(-u>) (25) 



with 



and 



B(u)B(-u) = uj 2 {uj 2 - F ) + {auj 2 - b 3 ) 2 (26) 



F 2 = lu 2 + 7^(1 — A) 
a =(l-j/A)n + 7 
b 3 =ttu 2 {l-A) 

Unlike the expression (fITjh the power spectrum (|25|) is represented by the ratio of two 
polynomials in u 2 , a quadratic one and a cubic one. Q 2 in the numerator may be 
neglected compared to u 2 around the resonance, if Q <C uoq. Simplifications in the 
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denominator will depend on the set of parameters A, y = 1 — x and 7 that are chosen, 
and to which Q/ouq must be compared. 

To link the present formalism to the previous one, it is appropriate to take 
AVL = Guq, because in the limit of high Q, small Q/ujq and for t^ujo = tt/2, the 



2 



M 0.5 - 



b£: 

: 




Figure 6. Comparison of analog protocol Equation (pj)) (the o points) and digital 
protocol Equation (|17j) with tdUio — n/2 (solid lines). The normalized power spectrum 
in Logarithmic scale in base 10 for two different settings of CI/uq: (a) = 0.2: (b) 
= 10~ 3 . Q = 10 5 ; x = Li n /L = 10 -2 . For each plot, settings from top to bottom: 
G = 0.06, G = 0.15, G = 0.4, G = 0.75. A in Equation ([25]) has been varied in all 
curves to implement the constraint AD, = Gloq. 
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definition Id (a;) = f^j-Ts(^) implies ~ ^p-I s (t — t d ). For the digital protocol 

we had Id(t) = GI s (t — td), and thus the identification. 

In Figure [6] we compare Equation ( jTTl) and Equation ( l25i) in the two cases of f2/a>o 
large and small, as we are now in the position to compare the two protocols. 

In Figure 6(a), Equation (1251) with Q/u = 0.2 is compared with Equation i fTTj) . 
with A = 0.3, A = 0.75, A = 2 and A = 3.75 (G = 0.06, G = 0.15, G = 0.4 and 
G = 0.75 respectively). Closer inspection reveals that with the analog protocol the 
shift of the resonance frequency, quantified by Au = io r — ujq, is equal to at the two 
extremes of the interval A e [0, 1], and is both non monotonic and negative inside the 
interval. A minimum of approximately Auj/ujq = —5 x 10~ 5 is attained in the vicinity 
of A = 0.5, before that quantity becomes monotonically increasing with A. In Figure 



6(b) we used Vl/luq = 10 for the analog protocol, again in comparison with Equation 



(DID, for the values A = 60, A = 150, A = 400 and A = 750 (i.e. G = 0.06, G = 0.15, 
G = 0.4 and G = 0.75, as before). In this case, the smallness of VL/ujq leads to an almost 
perfect overlap for all values of G. For this choice of fl/uio, the shift Au; appears to be 
consistently monotonic with growing A (and therefore G), in analogy with the digital 
protocol. 

As it transpires, the two systems are very similar for extremely high quality factors 
(here Q = 10 5 ), provided Q is sufficiently small. In fact, one could show that they 
overlap completely when we take Q/ljq — > (while keeping the ratio AQ/ujq constant) 
and we tune the gain factor by adjusting the amplification A to achieve the required 
gain G. This is somewhat surprising if we compare the expressions (TTTjl and ( 1251) . and 
it demonstrates the fact that an ideal filter has exactly the same effect on the current 
as that of a feedback with an ideal delay line. 

Finally, as we show in Figured the qualitative difference between the two feedback 
protocols starts to emerge when we consider low quality factors Q (for example equal 
to 1). From the power spectra derived from Equation ( 1251) one does not see the same 
features that can be seen for the digital protocol (Figure [5]), with the development of 
secondary peaks. With the analog protocol, no new characteristic oscillatory behavior 
emerges. The power spectra here are much more similar in shape to one another. 

Another feature that one is able to appreciate from Figure [7] is that for a given G 
a higher cooling efficiency is obtained by choosing A large and Q/uq small. 



6. A case study - the AURIGA detector 

The GW detector AURIGA is an example of the high Q case with the analog 
protocol. The detector can be modeled by three coupled low-loss resonators [H]: two 
mechanical ones (the bar and a plate of the capacitive transducer) and an LC electrical 
one [19]. Their dynamics is described by three normal modes at separate frequencies, 
and each mode is modeled as a RLC series electrical oscillator. Each mode is driven by 
Johnson noise and the resulting current is fed back to reduce its jitter below that set 
by the thermal bathjl]. The lowest frequency mode out of the 3 is well separated in 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Figure 7. Comparison of Equation (|25|) and Equation (fT7|) at low quality factor. 
The normalized power spectrum is shown in linear scale. Q = 1; x = Li n /L = 10 -2 ; 
G = 0.5. Solid line: phase-shifted feedback protocol Equation (|25|) . with = 10 -3 . 

X-dotted curve: Equation (f25|) with f^/t^o = 0.2. Dashed line: time-delayed feedback 
protocol Equation (|17p . with i^o = 7f/2. A in Equation (|25[) has been varied in all 
curves to implement the constraint AQ = Gwq. 



frequency from the other two and is the best approximation of a single oscillator [H [16]. 
Here fi/cu = 0.23, Q ~ 10 6 , x ~ 10~ 2 and 67 ~ 1.7 x 10~ 2 . Two different approaches 
have anticipated this scenario. One shows that, as in equilibrium, the power spectrum 
of the current is Lorentzian like, but with effective damping much higher (and effective 
temperature much lower) than in equilibrium [1]. The other approach verified that the 
probability distributions of the thermodynamic observables can be very well fitted by an 
equilibrium-like Langevin equation, with a modified temperature and modified damping 
governing the thermal driving [T6]. Given the high Q, the work in [16] approximated the 
analog protocol used in the experiment with a digital protocol in analyzing the stochastic 
time evolution of the system. Here we show that the approximation is valid also in the 



frequency domain as shown in Figure 6(a) and related discussion. At the same time 



the present approach highlights some nontrivial features pertaining to the effect of the 
feedback on the resonance frequency. 

7. Conclusions 



While it is customary to approach the theoretical analysis by working in the time 
domain with ideal delay lines, we have shown that in some situations there is a distinction 
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between obtaining an active feedback via this method and via filtering the frequencies. 
For high quality factors the difference is subtle, and can be nullified by appropriate 
choices of the settings. Conversely, there are irreducible differences for the responses to 
the feedback of harmonic oscillators with low quality factors. 

Our analysis improves our understanding of some nonequilibrium aspects of 
oscillators, currently of much interest in large deviations theory [20]. Lack of a complete 
thermodynamic description of the source and the feedback apparatus is still a problem in 
these contexts, since they are mostly concerned about dissipation, work and exchanged 
heat [21]. Using an instantaneous Langevin equation as a model, one typically studies 
the fluctuations of the injected power instead [221 [16], but its relation with the entropy 
production rate is subtle. Analysis of Langevin equations with explicit memory terms 
might offer insight and a new direction for approaching this problem. 

As a concluding remark, we wish to briefly discuss the problematic issues that one 
is likely to encounter if FDT-type relations were sought for in regard to the memory 
functions [15] . This attempt may well be suggested by the form of the Langevin equations 
fTTOj) and (|13p . essentially Langevin Equations with memory. Recall that FDT relations, 
in the presence of memory effects, assume causality conditions on the noise to be 
applicable, e.g. (N(t)I(t - t')) = for every t' > OpH [13]. However, from (HU) 
and (fT4j) . one sees that N(t), as well as N s (t), contain terms that are evidently causally 
related to I(t — t') and I s {t — t'), because they participated in the dynamical driving of 
the circuit at times earlier than t — t' (one just needs take ktj to be sufficiently large). 
In other words, any random thermal fluctuation of the bath exerts its influence on the 
measuring apparatus both instantaneously as well as at later times, via the recycling of 
the feedback. A similar analysis can be performed also with the second type of setting. 
The approach that we have used here seems crystalline and controllable, being based on 
the sole assumption on the nature of the self correlations of the thermal fluctuations. 
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